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ABSTRACT 

Dwarf, irregular and infrared-luminous starburst galaxies are all known to have 
"steep" luminosity functions with faint-end behavior roughly (j){L) oc L^^'^. This 
form is exactly what is expected if the luminosities of these objects fade with time as 
L oc t^^'"^, because the objects spend more time at low luminosities than high, even 
if they form with a wide range of initial masses. Models of young stellar populations 
show this fading behavior when the star formation has occured in a single, short, recent 
burst. Steep luminosity functions therefore do not require steep mass functions if the 
galaxies are powered by fading bursts. The local-galaxy Ha luminosity function — which 
is less steep than L~^'^ — is also well-fit by this mechanism, because ionizing photon 
flux fades much more quickly than broad-band optical luminosity. An age-luminosity 
relation and a wavelength-dependence of the luminosity function are both predicted. 
In the context of this mechanism, the slope of the luminosity function provides a 
constraint on the stellar initial mass function in the bursts. 

Subject headings: infrared: galaxies — galaxies: luminosity function, mass function 
— galaxies: starburst — galaxies: stellar content — stars: luminosity function, mass 
function 

1. Introduction 

While the normal field galaxy luminosity function (GLF) (j){L) (number density per unit 
luminosity) or <j){log L) (number density per unit log luminosity) is "flat" in optical bandpasses at 
the faint end, i.e., (t>{L) oc L^^ or 0(logL) = constant (Efstathiou, Ellis & Peterson 1988; Loveday 
et al 1992; Mobasher, Sharpies & Ellis, 1993; Marzke et al 1994a; Lin et al 1996, 1997; Gardner 
et al 1997; Ratcliffe et al 1997), many studies have found that objects in which the luminosity 
is thought to be dominated by young stars show a "steep" GLF, with roughly (/)(L) oc L~^-^ 
at the faint end. Parameterizing (l){L) oc L", the 60 fim. GLF from the IRAS Bright Galaxy 
Sample appears to show a = —1.8 at the faint end (Soifer et al 1987; although see Saunders et 
al 1990) despite the fact that these same objects lie in the flat part of the optical GLF. The 
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60 /xm luminosity is thought to originate in dust heated by the radiation from young stars at 
ultraviolet wavelengths where dusty galaxies are optically thick. Although the faint end of the 
optical GLF is flat, there may also be a small "upturn" at the very faintest end, around absolute 
magnitude Mb = —16 mag (Marzke et al 1994a; Driver & Phillips 1996; Loveday 1997), which is 
explained by a luminosity function with a = —1.8, among dwarf (i.e., low-luminosity) galaxies. 
In the case of the CfA survey, the upturn can be explained entirely by the luminosity function of 
the Sm-Im galaxies (identified on the basis of morphology) which show a = —1.87 it 0.2 (Marzke 
et al 1994b). Local galaxies spectrally classified as strongly or recently star- forming also show 
a steep luminosity function (Heyl et al 1997). A steep upturn at the faint end of the GLF is 
observed for dwarf galaxies in rich clusters with a ranging from —1.4 to —2.2 (Sandage, Binggeli 
k, Tammann 1985; Driver et al 1994b; Bernstein et al 1995; De Propris et al 1995; Lobo et al 
1996; Wilson et al 1997). A recent measurement of the luminosity function of dwarf galaxies or 
"knots" formed in the tidal tails of merging galaxies finds a = —1.75 it 0.27 in the i?-band for 
these objects (Hunsberger, Charlton &: Zaritsky 1996). Compact "super star clusters" observed in 
the vicinity of starburst galaxies or galaxy mergers and interpreted as the progenitors of globular 
clusters show a luminosity function consistent with a = —1.8 although the numbers are small 
(Lutz 1991; Holtzman et al 1992; Whitmore et al 1993; Conti &; Vacca 1994). Finally, a luminosity 
function with faint end behavior a —1.8 is often invoked as a natural explanation of faint galaxy 
counts and redshift distributions (Broadhurst, Ellis & Shanks, 1988; Ealcs 1993; Koo, Gronwall &; 
Bruzual 1993; Driver et al 1994a; Treyer & Silk 1994; Metcalfe et al 1995; Small et al 1995; Lilly 
et al 1995; Ellis et al 1996). In these studies, the steep GLF is largely required to account for the 
large numbers of faint blue galaxies, which are mainly irregulars (Glazebrook et al 1995a; Driver 
et al 1995; Abraham et al 1996) and thought to have luminosities dominated by young stars. 

The steep faint end of the GLF is usually attributed to a steep underlying galaxy mass function. 
A steep mass function at small halo mass is natural for cold and mixed dark matter models. 
In the Press &; Schechter (1974) formalism, on small mass- scales 0(M^) dM^ oc M^^^^""^/^^ dM^, 
where the post-recombination power spectrum of density fluctuations has P{k) oc /c", with n — > —3 
for adiabatic fluctuations on small scales. Thus (f>{Mh) dM^ oc Mj^'^dM^. This has been amply 
verifled by numerical simulation for both cold (Brainerd & Villumsen 1992) and mixed (Ma &: 
Bertschinger 1994) dark matter halos. However, the ejection of gas by early generations of stars 
in shallow potential wells implies that the mass converted to stars rises faster than linearly with 
M/j, so for identical stellar populations, the galaxy luminosity function should be shallower than 
the halo mass function (see Silk & Wyse 1993 for a review). Furthermore, in the IR- luminous 
galaxy sample, the large scatter in optical-IR colors (Soifcr et al 1987) and the lack of correlation 
between IR luminosity and galaxy mass inferred from rotation curves (Lehnert &i Heckman 1996) 
suggest that the starburst GLF is not strongly tied to the host galaxy mass function. 

In this Letter, we remark that there is a natural mechanism which ensures a steep GLF among 
young objects: Even if the galaxy mass function (where now by "galaxy mass" is meant "the 
mass of that part of galaxy's baryonic mass which is turned into stars") is flat, a GLF of roughly 
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the form (/){L) oc L~^-^ will be observed among any population of objects whose luminosities are 
dominated by light from short, recent bursts of star formation with a Salpeter-like initial mass 
function. This is because their luminosities decrease with time in such a way that they spend more 
time (and are therefore more numerous) at low luminosities than at high luminosities. This kind 
of mechanism underlies the theoretical explanation of the 60 /xm GLF by Scoville &; Soifer (1991) 
and a discussion of the super star cluster luminosity function by Meurer (1995). An important 
feature of this mechanism is that steep mass functions are not required to explain steep luminosity 
functions. 



2. Exposition 

The luminosity L of a main sequence star of mass M > IMq scales as L oc M^, where 
r] = 3.9 (Kippcnhahn & Wcigcrt 1990). So the main sequence lifetime r oc M/L oc M^'"^. Let an 
ensemble of stars be formed in a single burst with initial mass function (f>{M) dM oc M'^^"^) dM 
(where x = 1.35 is the Salpeter slope). The stages of stellar evolution beyond the main sequence 
have lifetimes and an integrated consumption of nuclear fuel less than or comparable to the main 
sequence values (Renzini & Buzzoni 1986). Thus for x < rj, the total luminosity of the evolving 
ensemble is dominated by stars just leaving the main sequence, with mass M = Mt^mit)- 

Ltot ^ C / (t){M)L{M) dM ~ M(j)(M)L(M) 

J\Mp, ri — X 



(1) 

Mturn(t) 



where C is a factor of order 2 (Renzini &: Buzzoni 1986) which accounts for the post-main sequence 
fuel consumption. Since Mturn(0 is defined hy t = T(Mturn)5 we have for x < r), 

LTOTO^t^, C=^, (2) 
rj — 1 

so in this approximation ~ —0.9 for a Salpeter IMF. The luminosity in a particular wavelength 
band X will have a somewhat different dependence on age than the bolometric luminosity, 
depending on bolometric correction, but since stars consume about as much fuel in the (red) 
post-main seqence as during the (blue) main sequence, the dependence in optical and near IR 
bands will not be vastly different from that of the bolometric luminosity. 

If in wavelength band X an object has luminosity Lx which varies with time t since its birth 

as 

Lx oc (3) 

and the object is observed at a random times after its birth, the probability that it is measured 
to have luminosity Lx is proportional to {dLx / dt)~^ oc L^^^^'' ^ . Therefore, a population of 
identical objects of this type born at a fairly uniform rate over a time interval will appear to have 
a luminosity function with faint end slope ax = (1/Cx) — 1 despite being intrinsically identical. 
The luminosity function is produced by a spread in ages rather than a distribution of intrinsic 
properties. To make a luminosity function with a —1.8, C, ~ —1.3 is required. 
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In any realistic scenario, the objects will have a distribution of intrinsic sizes or starburst 
masses. Consider objects with luminosities Lx = M Kx{t) where M is the mass of stars created 
in the starburst, the objects are selected from a distribution 0(-/Vf) of starburst masses, and 
Ax(t) oc t^^ is the light per unit mass in wavelength band X as a function of time. The faint 
end of the luminosity function will still show 4'{L) cx L^^^^'>^'^ as long as 4){M) is flatter than 
<p{M) oc M^^/'^')^^ (i-e., the exponent is less negative than (1/C) — 1) at the small-mass end and 
steeper than oc M^-^^''')^^ at the large-mass end. Both conditions hold if the "flat" Schechter 

function observed for the local GLF represents the intrinsic mass distribution for galaxies. 

With models much more detailed than the scaling argument given at the beginning of this 
section, Leitherer & Heckman (1995) and Bruzual &: Chariot (1993) predict the fading with time of 
young stellar populations in a number of photographic bands and for several models of the stellar 
initial mass function (IMF) and metallicity, to which power laws have been fit. The results are 
given in Table ||. The exponents in the ultraviolet and optical, for Salpeter IMF (slope x = 1.35) 
are in the range —1.4 < C < —0.9, which lead to luminosity functions with —2.1 < a < —1.7. A 
Miller-Scalo IMF (slope x = 2.3) produces slower fading and therefore steeper (more negative a) 
luminosity functions. Thus if the luminosity function of faint galaxies is indeed produced by fading 
of young populations, the slope of the luminosity function can be used to constrain the stellar 
IMF. The different luminosity function slopes observed in different clusters (Binggeli, Sandage &: 
Tammann 1985; Driver et al 1994b; Bernstein et al 1995; De Propris et al 1995; Lobo et al 1996; 
Wilson et al 1997) could reflect an environmental dependence in the stellar IMF, although the 
statistics are small enough that the various studies appear to be at least marginally consistent 
with one another. 

In addition, the evolution of the rate Fj of production of ionizing photons (i.e., photons 
shortward of 912 A) is predicted (Leitherer & Heckman, 1995), with roughly Fj oc t^^'^. If the 
strength of Ha emission is proportional to Fj, this fading behavior naturally produces an Ha 
luminosity function with faint-end slope a = —1.2, consistent with a recent measurement of 
a = -1.3 ± 0.2 for the local Ha GLF (Gallego et al 1995). 

As an illustration. Fig. 1 shows: (a) fading laws of the form 

Aboi(t) = for t < to 

Ao (^)'^ fort>to (4) 

with = —1.4 and —0.7 which serve as approximations to Fig. 7 of Leitherer & Heckman (1995); 
(b) distributions of starburst masses (/>(logM) (number per decade in mass) of the Schechter 
(1976) form 

0(log M) a M 0(M) oc i^—j e"^^/*^* (5) 

where /3 sets the slope of the mass function at the low-mass end and M* (fixed at 10^ Mq for 
definiteness) sets the high-mass cutoff; (c) the resulting luminosity functions (j){logL) oc L4>{L) 
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(number per decade in luminosity) under the assumption that the populations are observed at 
random ages as they fade; and (d) the age distributions (/>(logt) (number per decade in age) at 
two different luminosity levels (10^^'^ and 10^^'^ Lq in this case) . This shows that the luminosity 
function has the predicted power-law form even if there is a broad distribution of starburst masses. 
It also shows that the distribution of ages can be broad but does indeed get older as the objects 
get fainter. The age function is broadest in the cases in which the mass function is almost as steep 
as the luminosity function. 

As an illustration, the upturn in the local GLF at absolute magnitude Mb ~ — 16 
{Hq = 100 kms^^Mpc"^) found by Marzke et al (1994a) can be analyzed in terms of the fading 
mechanism. The GLF has amplitude 0(logL) ~ 0.09 Mpc~'^ at absolute magnitude Mb = —16 or 
luminosity Lb = 3.7 x 10^ Lq and slope a ~ —2 faintward. To order of magnitude, by dimensional 
analysis, 

dAB 



(/.(log Lb) = (In 10) Lb = (In 10) Lb '^^ 



dt 



(6) 



dVdLB dV dt Mtyp 

where dN/ {dV dt) is the birth rate density (number per unit time per unit volume) , -Mtyp is 
the typical starburst mass, and dh^B/dt is the rate of change of the light-to-mass ratio in the B 
band, evaluated at the light-to-mass ratio = Lb /Mtyp at which starbursts of mass Mtyp have 
luminosity Lb- At this value. 



dAB CsAo ftX'^B-i Cb^o f L 



1-(1/Cs) 



dt to \to) to l^MtypAoj 

where Aq = 5.8 x Lq/Mq, to = 3 x 10^ yr and (^b = —1.0 are the appropriate values for the 
fading-law parameters in equation (^) for the B band (Leitherer & Heckman 1995). Plugging in, a 
relationship between the birth rate density and the typical mass is derived 

dN _\CB\H^ogLB) fMtypAo\y^B_ 3 / Mtyp 



N 1 . ; - 8.3 X 10'^ ] Mpc-^Gyr"^ (8) 

dVdt (In 10) to y Lb J \W^MqJ ^ ^ ^ ' 

for Ho = 100 kms^^Mpc^^. Note that the model is not unique because neither the birth 
rate density nor the typical mass can be determined from the above information alone, only a 
relationship between the two. 

Of course this fading mechanism requires a fairly uniform birth rate. In the extreme, if in 
some population all starbursts occured simultaneously at some moment in the past, no spread of 
ages would be observed, and the luminosities would be proportional only to the intrinsic sizes of 
the starbursts. The GLF would have nothing to do with the fading exponent C,x at all. If the birth 
rate is varying slowly with time, it will affect the final luminosity function shape. For instance, if 
the birth rate is decreasing with time, there will be more old galaxies than young and therefore 
more low-luminosity systems than would be predicted for constant birth rate and a steeper GLF 
will be observed (i.e., a will be more negative than [l/Qx) — 1)- Again, different slopes among, 
say, different galaxy cluster GLFs could indicate different dependences of starburst birth rate on 
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time. However, over the small dynamic ranges (and hence age ranges) over which such GLFs have 
usually been measured, the assumption of relatively constant birth rate is not unreasonable. 

If individual galaxies undergo multiple starbursts, this fading mechanism will still ensure a 
steep GLF, as long as the sources are fading according to Ax{t) oc t^-^. For a multiple-burst source 
this will be true at young ages when the recently formed stars still dominate the light. However, 
at late times the underlying population of old stars from previous bursts will eventually become 
a significant, non-fading contribution to the luminosity. So multiple-burst sources will show the 
steep GLF, but perhaps over a somewhat restricted range of age and therefore luminosity. Of 
course because young stellar populations are so much brighter per unit mass than old, this will not 
be a big restriction unless the sources burst very frequently (in which case the basic assumption of 
the model breaks down). 

3. Prediction 

This fading model for the faint end of the luminosity function makes three important 
predictions. The first prediction is that the bright members of the populations which show these 
steep luminosity functions should be mostly young. The distribution of ages and in particular 
the modal age tm{Lx) at each luminosity Lx in band X should obey the law tm{Lx) oc L}^^^ . 
A correlation between dust temperature and luminosity is observed in IRAS galaxies (Miley, 
Neugebauer &: Soifer 1985; Rieke & Lcbofsky 1986; Soifer et al 1987). In models in which the 
far-infrared emission is from dust heated by young stars, this relation can be interpreted as an 
age-luminosity relation, although it is also naturally produced if the dust is heated by other 
mechanisms (e.g., nuclear activity or continuous star formation). Spectra of IRAS galaxies in the 
0(L) oc L-i-S part of the 60 /xm luminosity function could be used to confirm this interpretation 
of the temperature-luminosity relation. In the optical an age-luminosity relation for a starburst 
can be converted into a color-luminosity relation; the expectation is that fainter galaxies will be 
redder, in possible contradiction to the evidence from some samples of cxtragalactic HII regions 
(Telles 1995) and faint galaxy counts (e.g., Koo &: Kron 1992). Again, spectra of nearby dwarf, 
irregular or infrared-luminous galaxies could in principle be used to determine a quantitative 
age-luminosity relation. Because the luminosities are determined more by age than by intrinsic 
properties, there ought to be little correlation between the luminosities and dynamical masses of 
the galaxies, as measured from rotation curves or velocity dispersions. Of course the width of 
the age distribution at a given luminosity depends on the underlying starburst mass distribution 
discussed above. If the starburst mass distribution 0(m) is a Schechter function with faint end 
slope (5 the width of the age distribution increases as /3 — a = (1/C) — 1; see Fig. l('dj. (Recall 
that if /3 < q;, the overall luminosity function will have faint-end slope (3 rather than a.) In 
principle this effect could be used to constrain the distribution of starburst masses. 

The luminosity function shape depends on the objects fading out of the sample with time, so 
the second prediction is that this mechanism can only explain the slope of the luminosity function 
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over as many magnitudes as the fading obeys the power law. This range appears to be at least 
7.5 mag (factor of 10^) for the visible-band luminosities (Bruzual & Chariot, 1993) — unless the 
decay is interrupted by a subsequent burst of star formation, of course — but luminosity functions 
with steep faint-end slopes over a significantly greater range cannot be explained entirely by 
this mechanism. The largest dynamic range for which the slope has been measured so far is 
in the IRAS sample which appears to show slope a = —1.8 over a factor of roughly 300 in 
60 /im luminosity (Soifer et al 1987). The luminosity functions of cluster galaxies will have to be 
measured to levels more than three magnitudes fainter than current observational limits in order 
to exceed the allowed range. 

The third prediction comes from the different slopes of A(t) in different wavelength bands. 
These translate into different predicted luminosity function slopes a. In general, the luminosity 
function is expected to become steeper (a more negative) as the objects are observed in 
longer- wavelength bands. This prediction is related to the color-luminosity relation predicted 
above and is easy to test with multi-band imaging of dwarf galaxy populations. There is already 
some evidence in contradiction to this prediction for galaxy cluster GLFs (Wilson et al 1997), and 
faint galaxy counts, which may be naturally explained with steep luminosity functions, become 
flatter, not steeper, with increasing wavelength (e.g., Koo &: Kron 1992). Of course the faint 
galaxy counts are difficult to interpret, since the observed objects are effectively integrated over a 
range of redshifts (and hence rest wavelengths) and there is the possibility that multiple distinct 
populations play comparable roles. It is encouraging that the Ha luminosity function (Gallego 
et al 1995) is well-fit by this mechanism; under this interpretation the Ha luminosity becomes 
more a measure of time since most recent starburst than a measure of star-formation rate, as is 
conventionally assumed. Conclusions about the volume averaged star-formation rate of the local 
Universe (Gallego et al 1995) are not strongly affected by this change in interpretation, however, 
because such conclusions only depend on total emission of ionizing photons, not whether or not 
the emission is simultaneous with or subsequent to the star formation. 

4. Application 

Finally, we emphasize that if the fading of short, recent bursts of star formation does indeed 
explain the faint end slope of galaxy luminosity functions, the measured slope becomes a strong 
constraint on the stellar IMF in these objects. Since other techniques (e.g. population synthesis 
models) often do not constrain the IMF uniquely in individual objects (e.g., Santos et al 1995), 
the slope of the faint end of the galaxy luminosity function may become an extremely useful 
diagnostic for star and galaxy formation models. 

We benefitted from helpful comments from Lee Armus, Roger Blandford, Tim Heckman, 
Jeremy Heyl, Nick Scoville, Tom Soifer and an anonymous referee. Glaus Leitherer, Tim Heckman 
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Table 1. Exponents of Luminosity Decay' 
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= 1.35, Z 
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-1.2 


-1.0 


-0.9 


-0.7 


Ref 2": 
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= 1.35, Z 


= O.lZo 


-1.4 


-1.2 


-1.0 


-1.0 


-0.5 


Ref 2": 
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= 2.3, Z = 


-- Q.IZq 


-0.7 


-0.5 


-0.4 


-0.3 


-0.0 



^See equation (g). 

''Ref 1: Bruzual & Chariot 1993. Luminosities fitted over the age 
range lO^y < t < 10^'' '*y. IMF slope x (equation |l|) relevant only for 
stars with masses greater than 1 Mq. 

"^Ref 2: Leitherer & Heckman 1995. Luminosities fitted over the age 
range 10^'^y < t < 10^'^y. IMF slope x (equation 0) relevant only for 
stars with masses greater than 2Mq. 
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Fig. 1. — An illustration of the fading mechanism for the steep GLF presented in this Letter. Shown 

are (a) the bolometric fading laws Aboi(t) (bolometric luminosity per unit mass in solar units) used, 
crude fits to Fig. 7 of Leitherer & Heckman (1995); (b) the mass functions (p{\ogM) oc M (j){M) 
(number per decade in mass), really the distributions of starburst masses, arbitrarily choosing 
M* = 10^ M©; (c) the bolometric luminosity functions 0(logLboi) oc Lboi0(Lboi) (number per 
decade in luminosity) produced by the fading and mass functions; and (d) the age function 
^{iogt) oc t(/)(t) (number per decade in age) for objects at luminosity Lboi = 10^^'^ (lower 
curves) and Lboi = lO^*^'^ Lq (upper curves). Three models are shown, one with C, = —1.4 and 
/? = -1.0 (solid line), C = -1-4 and (3 = -1.75 (dotted), and C = -0.7 and (3 = -1.0 (dashed). See 
text for symbol definitions. Mass, luminosity and age functions all have arbitrary normalizations; 
it is not meaningful to compare relative heights of two different curves on the same plot. 



